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Black holes, an extreme consequence of the mathematics of General Rel- 
ativity, have long been suspected of being the prime movers of quasars, 
which emit more energy than any other objects in the Universe. Re- 
cent evidence indicates that supermassive black holes, which are probably 
quasar remnants, reside at the centers of most galaxies. As our knowledge 
of the demographics of these relics of a violent earlier Universe improve, 
we see tantalizing clues that they participated intimately in the formation 
of galaxies and have strongly influenced their present— day structure. 

CO 



Black holes are a prediction of Einstein's theory of gravity, foreshadowed by the work 
of Michell and later Laplace in the late 18th century. K. Schwarzschild discovered 



the simplest kind of black hole in the first solution of Einstein's equations of General 



Relativity, and Oppenheimer was among the first to consider the possibility that black 
holes might actually form in nature. The subject gained life in the 1960s and 70s, 
when supermassive black holes were implicated as the powerhouses for quasars and 
stellar-mass black holes were touted as the engines for many galactic X-ray sources. 
In the last decade, we have progressed from seeking supermassive black holes in only 
the most energetic astrophysical contexts, to suspecting that they may be routinely 
present at the centers of galaxies |1], |2|, |3|] . 

The defining property of a black hole is its event horizon. Since the horizon itself 
is invisible, we must often settle for evidence of mass without light. All dynamical 
techniques for finding supermassive black holes at the centers of galaxies rely on a 
determination of mass enclosed within a radius r from the velocity v of test particles; 
in Newtonian physics, this mass is M r = a v 2 r / G. Determining a requires a detailed 
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dynamical analysis, but it is often of order 1. In cases where there is extra mass above 
that associated with starlight, we refer to the object as a "massive dark object" (or 
MDO). In most of the cases discussed in this paper, it is likely that the MDO is a 
supermassive black hole (MBH), but in only a few cases have plausible alternatives 
to a black hole been ruled out. These are important as they establish the reality of 
MBH and justify the interpretation of less compelling objects as MBHs. 

1 Black holes as the energy sources of quasars 

Black holes are thought to exist in two mass ranges. Small ones of ~ 10M Q are 
the evolutionary end points of some massive stars. This paper discusses the much 
more massive ones that might power quasars and their weaker kin, active galactic 
nuclei (AGN). Quasars produce luminosities of L ~ 10 46 ergs/sec (~ 10 12 L Q ). Where 
they power double-lobed radio sources, the minimum energy stored in the lobes is 
E ~ 1Q60-64 er g S _ ^he m ass equivalent of this energy is M = E/c 2 ~ 1O 6 ~ 1O M , 
and the horizon scale associated with that mass is Rs = GM/c 2 ~ 10 11-15 cm. 
Although most quasars do not vary much at visual wavelengths, a few objects change 
their luminosity in minutes at high energies @, Since an object cannot causally 
vary faster than the light-travel time across it, such objects must be smaller than 
R ~ cr ~ 10 13 cm. Although relativistic corrections can alter this limit somewhat in 
either direction via Doppler boosting or gravitational redshift, there is no escaping 
the conclusion that many quasars are prodigiously luminous yet tiny, outshining a 
galaxy in a volume smaller than the solar system. 

The small size together with the enormous energy output of quasars mandates black 
hole accretion as the energy source. Most investigators believe that quasars and AGNs 
are MBHs accreting mass from their environment, nearly always at the center of a 
galaxy || |7], || . Black holes of mass > 1O 7 M must normally lie at the center because 
dynamical friction drags them to the bottom of the potential well. This location is 
now clearly established for low redshift (z <0.3) quasars ||. The connection between 
MBHs and quasars was first made by Zeldovich [[3IJ and Salpeter pl| . Lynden-Bell 
fT2|| sharpened the argument by computing the ratio of gravitational energy to nuclear 
energy 



where Rs is the Schwarzschild radius of a black hole of mass M, R is the size of the 
quasar, and e g and e n are gravitational and nuclear energy conversion efficiencies; the 
last equality follows from the typical astrophysical thermonuclear efficiency of ~ 1% 
and the size scale from variability noted above. 

Because quasars were populous in the youthful Universe, but have mostly died out, 
the Universe should be populated with relic black holes whose average mass density 
p u matches or exceeds the mass-equivalent of the energy density u emitted by them 
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13] . The integrated comoving energy density in quasar light (as emitted) is 

dt 



u 
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/ / L§(L\z)dL — dz = 1.3 x 10~ 15 erg cm~ 3 , (2) 
Jo Jo dz 



where $ is the comoving density of quasars of luminosity L, and t is cosmic time. The 
corresponding present-day mass density for a radiative efficiency e is p u = u/[ec 2 } = 
2 x 10 5 ^ M Mpc~ 3 . This density can be compared to the luminous density in 
galaxies, j = 1.1 x lO 8 L Mpc~ 3 0, to obt ain the ratio of the mass in relic MBHs 
to the light of galaxies: 



T =7 =L8X1 ° ItA^J- (3) 



2 Dynamical evidence for massive black holes 
2.1 First steps 

The first dynamical evidence for black holes in galactic centers was the 1978 measure- 
ment ]T3 of a rising central velocity dispersion, reaching ~ 400 km s _1 , in the giant 
elliptical galaxy M87. This object is a prime site to prospect for an MBH by virtue 
of its AGN features — nonthermal radio emission, broad nuclear emission lines, and 
a "jet" of collimated relativistic particles being ejected from the nucleus. Isotropic 
models of the stellar kinematics, when combined with photometry, implied an MDO 
of 5 x 1O 9 M . The result was criticized because the data were also matched by a 
model with radially anisotropic stellar orbits and no black hole. Thus, the importance 
of understanding the stellar orbital structure of the centers of galaxies was obvious 
at the very beginning, and this subject has developed in parallel with the search for 
MBHs. 

More convincing evidence was found in the 1980s for MDOs in M31 and M32 (the 
Andromeda galaxy and its satellite), which are nearby and hence observable at high 
spatial resolution |L6], [nj. Rapid rotation near their centers reduces the danger of 
confusing a central mass with radial orbits, and Schwarzschild's method of loading 
orbits in a galaxy potential ]TB, [ll| was used to eliminate unphysical models. Modern 



methods use Schwarzschild's method to fit the entire line-of-sight velocity distribu- 
tion for axisymmetric models. 

The stellar velocity work on M87 was largely vindicated two decades later by the 
Hubble Space Telescope (HST), which revealed a small gas disk at the center |20], pT| . 
The gas is plausibly in circular motion, so the MDO mass estimate is straightforward. 



These and later data provide strong evidence for an MDO of 3 x 10 M [20, 22 



value similar to but slightly smaller than the one derived in [23 



3 



2.2 Two remarkable examples 



The work described above revealed strong examples of MDOs but no iron-clad evi- 
dence for MBHs. This gap has now been partly closed in two remarkable objects. The 
mild AGN NGC 4258 was shown to possess a tiny annular gaseous disk near the 
nucleus, populated by water masers whose Doppler velocities can be observed with 
exquisite precision. The rotation curve is Keplerian to high accuracy over the annulus 
width (0.13 — 0.26 pc). The very small velocity residuals of "zl% inspire confidence 
in the derived mass of 3.6 x 10 7 M© The extraordinarily high implied density 

of > 10 9 M pc~ 3 (10 12 if one takes the limits on departures from Keplerian motion 
as a constraint on the concentration of the mass) permits the use of astrophysical 
arguments to rule out most other explanations for the dark mass (see below). This 
is a firm link from MDOs to MBHs. 

The center of our Galaxy holds the second confirmed MBH. Near-infrared observa- 



tions detect proper motions of stars in orbit about the galactic center [gt], [27|] and 
indicate a rising stellar velocity dispersion down to distances of 0.01 pc. For the first 
time stars are being observed to orbit an MDO, year by year, with impressive accu- 
racy that will steadily improve. The density of > 10 12 M & Mpc -3 within the resolved 
region is again extraordinarily high, ruling out most alternatives to an MBH. 



In the Galaxy and NGC 4258, the MDOs are almost surely MBHs [2l| rather than 
clusters of smaller masses. A cluster of radius r and total mass M of N self-gravitating 
point masses will collapse or evaporate on a timescale of a few hundred two-body 
relaxation times, t r = O.UN (r 3 /GM) 1/2 [h^O^iV)]- 1 (some important caveats were 



noted in PS|). At a chosen nonrelativistic density, the lifetime of the cluster can 
be made longer than the age of the Universe by making the point masses sufficiently 
light (and therefore numerous). In NGC 4258 and the Galaxy, this constraint requires 
constituents with masses <O.1M . Brown dwarfs and white dwarfs with this mass 
have large radii and would rapidly collide and merge. Thus the remaining candidate 
components for cluster models are low-mass black holes or noninteracting elementary 
particles. However, there is no known way to make O.1M black holes, and, since 
non-interacting elementary particles do not radiate energy efficiently they do not 
settle into a compact configuration. The cases for black holes in both NGC 4258 and 
the Galaxy therefore seem very strong. The argument for all of the other MDOs is 
weaker, but, by analogy and by virtue of the overwhelming circumstantial evidence 
for MBHs in AGNs they are now the preferred explanation. 



2.3 Potential new tools 

The techniques above are difficult to apply to galaxies containing AGNs. The bright 
nucleus renders the stellar absorption lines close to the center nearly invisible, and 
in many cases the nebular emission lines are influenced by nongravitational forces. 
Nevertheless, in AGNs that emit broad emission lines originating from gas near the 
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central engine (Aw > 10, 000 km s _1 ), one can attempt to estimate a mass from the 
average velocity of the gas and the radius of the emitting region. The velocity comes 
from the widths of the lines, but the radius is harder to measure. It can be estimated 
either from photoionization models of the gas or by "reverberation mapping." In 
the latter method, the radius is inferred from the time interval (due to light travel) 
between fluctuations in the continuum radiation and the changes these induce in the 
emission lines |29| . 



Masses obtained from reverberation mapping for nearby Seyfert nuclei range from 
10 7 to 10 8 M Q [30|, roughly consistent with (but somewhat smaller than) the dynam- 



ical results of the previous section. However, the identifications of line width with 
orbital velocity and time delay with radius are problematic given the absence of any 
correlation between line width and radius of the form v oc r -1 / 2 within the same 
object. Evidently some essential component of the model is still missing. A proper 
understanding of this technique would give us a powerful tool for more luminous and 
distant objects . 



Recent advances in X-ray astronomy have furnished dramatic new evidence for MBHs 
in AGNs. It had been known for some time that the X-ray spectra of many AGNs 
show an iron Ka emission line at a rest energy of 6.4 keV, thought to arise from X-ray 
fluorescence of cold, neutral material in an accretion disk. Until recently, the available 
spectral resolution was insufficient to test the predicted line profile, but the Advanced 
Satellite for Cosmology and Astrophysics provided the much-awaited breakthrough 
in the Seyfert 1 galaxy MCG-6-30-15 P^f . The Ka line exhibits relativistic Doppler 
motions of nearly 100,000 km s -1 , as well as an asymmetric red wing consistent 
with gravitational redshift. The best-fitting disk has an inner radius of only a few 
Schwarzschild radii. The Fe Ka line profile has now been seen in many objects [[33 



and data of better quality may eventually even allow measurement of the spin of the 
black hole M,M. 



3 The demographics of supermassive central black 
holes 

We now turn to the questions of the number of MBHs in the Universe. Figure 
1 illustrates the relationship between black hole mass and host spheroid luminosity 
from the data in Table 1 (labeled points). The labeled MDO masses seem to correlate 
with spheroid luminosity (solid line); the upper limits are also consistent with this 
relation 0, |36|| . However, the number of points is small, and, further progress requires 
more objects. Few objects at the present time have been studied with the detailed 
spatial resolution and/or modeling of the labeled points. 
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Mass estimates of the candidate MBHs in galaxies with dynamical information plotted 
against the bulge luminosity of their host galaxy. The labeled points are the results 
of painstaking observation and detailed modelling. The symbols indicate the how M, 
was derived: kinematics of gas — triangles; dynamics of stars — filled circles; masers 
— diamonds; or two-integral modelling using ground-based stellar kinematics — small 
squares. Arrows indicate upper limits on M.. The solid line is a model with M. = 
0.005M bu i ge and . M bu i ge = 5(L bu i ge /10 9 L Q ) 1 - 2 . The distribution of M, is roughly 
Gaussian in log(M./M bu ige) with mean —2.27 (M,/M bu i ge = 0.005) and standard 
deviation 0.5. The dashed line is the quasar light prediction of eqn 3 apportioned 
according to the bulge mass: M, = 2 x 10 7 (L bu i ge /5 x 10 9 L Q ) 1 - 2 . 
The small offset from the observed black- ho le/bulge-mass relation indicates that the 
present integrated density in MBHs is broadly consistent with the integrated lumi- 
nosity produced by AGNs over the life of the Universe. This offset may reflect a 
radiative efficiency of average quasar accretion less than 0.10. 
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DYNAMICALLY IDENTIFIED MDOs 



Galaxy 


Name 


Type 


Distance Q 


M B 


M. 


Reference 








Mpc 


Bulge 


M sun 





Stellar Dynamics 





Galaxy 


Sbc 


0.0085 


-17.70 


NGC0221 


M32 


E2 


0.7 


-15.51 


NGC0224 


M31 


Sb 


0.7 


-18.82 


NGC3115 




SO 


8.4 


-19.90 


NGC3377 




E5 


9.9 


-18.80 


NGC3379 


M105 


El 


9.9 


-19.79 


NGC4342 


IC3256 


SO 


15.3 


-17.04 


NGC4486B 




El 


15.3 


-16.65 


NGC4594 


M104 


Sa 


9.2 


-20.88 








Gas Dynamics 




NGC4374 


M84 


El 


15.3 


-20.96 


NGC4486 


M87 


EO 


15.3 


-21.42 


NGC4261 




E2 


27.4 


-20.82 


NGC7052 




E4 


58.7 


-21.31 



2.8 x 10 6 
3.4 x 10 6 
3.0 x 10 7 
2.0 x 10 9 
1.8 x 10 8 
6.7 x 10 7 
3.0 x 10 8 
5.7 x 10 8 
1.0 x 10 9 



sec 



mm 



0, 0, |74], [75 
[0,0 



©0 

H 





n n 



1.4 x io 9 ga 

3.3 x 10 9 

4.5 x 10 8 
3.3 x 10 8 || 



2C, 21 



8 J 



Maser Dynamics 
NGC1068 M77 Sb 15. -18.82 

NGC4258 M106 Sbc 7.5 -17.28 

NGC4945 Scd 3.7 -15.14 



i.o x io 7 m 

4.2 x 10 7 
1.4 x 10 6 



NGC0205 
NGC0598 
NGC4395 



M33 



Spheroidal 
Scd 
Sm 



Upper Limits 

0.72 -9.02 
0.795 -10.21 
-7.27 



< 9. x io 4 m 

< 5. x IO 4 [|U 

< 8. x IO 4 



"All quantities in this paper are computed for a Friedman-Robertson- Walker 
Universe with = 1 and Hq = 80 km s~ Mpc~ . Distances to nearby MBHs come 
from many sources, but are always rescaled to this Hubble constant. 

At the risk of greater uncertainty, more galaxies can be included by combining ground- 
based stellar kinematics with HST central light profiles. A simple modeling recipe 
based on two-integral axisymmetric models has been used on such data for a further 
25 E and SO galaxies [[57]. This procedure assumes that the phase-space density is only 
a function of the energy and one component of angular momentum. MDO masses from 
this technique can be checked against galaxies with HST spectroscopy. The results for 
5 low-mass galaxies with steep inner light profiles show good agreement, but the 2- 
integral method may overestimate masses by a factor of a few for massive ellipticals 
like M87, which rotate slowly and have shallow central light profiles. A Bayesian 
analysis of this sample indicates that MDOs are in fact very common features of 



normal, bright galaxy centers [37] . 



All known MDOs with measured masses so far are in galaxies with identifiable 
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spheroidal components, suggesting that black hole formation is exclusively linked to 
spheroid formation. However, several well-imaged low z quasars do not appear to be 
associated with spheroids ||. 

It is clearly important to survey more late-type spirals without bulges. So far, we 
have little dynamical evidence on the centers of such galaxies, but AGN activity might 
perhaps be used as a proxy for a black hole. Seyferts are generally not found in late- 
type spirals f38|, but a single dwarf Seyfert nucleus (out of hundreds surveyed) has 



been discovered in NGC 4395, a nearby bulgeless Sd IV galaxy The bolometric Seyfert 
luminosity of this nucleus is 1.4 x 10 40 ergs/sec |39|], and the Eddington BH mass is 



only 11OM , small enough to have been produced by stellar evolution. 

There is no evidence for MDOs in low-surface brightness galaxies, although existing 



studies do not set compelling limits [40 



Our view of these observational results, largely developed over the past decade, is 
as follows: (i) MBHs are a normal feature of the central regions of bright galax- 
ies, particularly those with spheroids; (ii) their masses scale in rough proportion to 
host-galaxy spheroid mass; and (iii) the total mass density in black holes is broadly 
consistent with the mass-equivalent energy density in the quasar light background. 
We therefore believe that the black hole fossils of the quasar era have been found. 



4 Co— evolution of Galaxies and Black Holes 



4.1 The Era of Quasars 

The improving statistics on local MBHs can be compared to their properties and 
distribution during the quasar era at z ~ 3. Were today's MBHs already fully formed 
by that time, or was the average MBH smaller then, having grown by later accretion 
or mergers to form the present population? The evidence is not conclusive but seems 
to favor some growth. 

The epoch of maximal quasar activity in the Universe peaked at the same time, or 
slightly before, the epoch of maximal star formation, and MBHs must have formed 
before this time to be available to power quasars. Figure 2 illustrates this point by 
plotting the history of the rate of observed star formation in the Universe BTI, [43| 



together with the density of luminous quasars (those with L > 3 x 10 46 ergs/sec, p4H ). 
The rise in starbirth is tracked closely by the rise of luminous quasars. However 
the bright quasars reach their peak at z>2 (t < 1.6 x 10 9 yr), and then proceed 
to die off nearly 10 9 yr before the peak in star formation, which occurs at z ~ 1.2 
(t = 2.6 x 10 9 yr) f4"5fl . The application of extinction corrections to the cosmic star 
formation rate before z ~ 2 is a controversial subject that may evolve rapidly as 
better IR data become available. 

This chronology favors models in which the black hole forms before, or in close asso- 
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Figure 2. — 

A comparison of the density of very bright quasars in the Universe, with the density 
of star formation, as a function of redshift and cosmic time. The solid line and the 
filled circles are the comoving number density of quasars (uq see [12]) in units of 
Mpc -3 , while the open circles represent estimates of the comoving star formation rate 
(SFR). The SFR can be read off the left axis, and n Q should be read off the right axis. 
The arrows on the two highest redshift SFR points indicate an estimate of plausible 
correction for extinction. The peak of quasar activity in the Universe appears to 
predate most of the star formation. The times are derived from the redshift assuming 
n = 1,H Q = 80 km s~7 M P c - 
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ciation with, the densest parts of galaxies (for example [fSJ), as opposed to models 
in which the galaxy forms first and later spawns the quasar. For this reason we as- 
sociate the birth of quasars with spheroid formation, a process also identified with 
dense regions that collapse early |47|]. The collapse of these regions would predate 
the average rise of star formation illustrated in Figure 2. 

The decline of bright quasars at z < 2 seems likely to reflect a loss of fuel. Galaxy 
mergers, an effective gas transport process, are less frequent as time passes and involve 
a lower mean density and slower dynamical clock. There is also less gas in galaxies 
overall (especially spheroids). Radial transport instead may be choked by the onset 
of chaos, which would make the spheroid axisymmetric [081 or through the creation 



of an inner Lindblad resonance |49|j. Any or all of these effects could play a role in 



the reduction of fuel for the black hole. 

A model in which an MBH forms in every collapsing spheroid invites comparison 
between the densities of luminous quasars and galaxies at z ~ 2, and of fossil MBHs 
today. To make this comparison, we must identify a quasar of specified luminosity 
with its fossil MBH mass. We do this under the assumption that the brightest quasars 
are Eddington-limited, that is, that their luminosity is so great that radiation pressure 
on nearby electrons balances the gravitational force on associated protons. In this 
situation, the "Eddington luminosity" is 

L E = 1.3 x 10 46 M 8 ergs/sec, (4) 

where M 8 is the mass in units of 1O 8 M . 

The comoving density of luminous quasars with L > 6 x 10 45 ergs/sec reached its 
peak value of 10~ 6 Mpc~ 3 at z ~ 3 [44]. We assume for the moment that the MBHs 



underlying these bright quasars have not grown substantially since then. To estimate 
the corresponding mass range, we first correct the absolute magnitude limit in the 



QSO surveys (Mg < —26.0, |44], |50fl ) downward to our H and upward by a bolo- 
metric correction factor of 10 |)T| to get a lower-limit luminosity of 6 x 10 46 ergs/sec. 
This luminosity implies a MBH mass of M ~ 4 x 1O 8 M , corresponding to a bulge 
luminosity of current hosts (from Figure 2) of ~ 10 10 L & . Correcting for the total-to- 
bulge luminosity ratio of 3 []5l|, and using the luminosity function of bright galaxies 
TM, we identify the bright quasars of z>2 with the half of modern galaxies with 



Mb < —20.7 that have bulges |53| , we find a comoving density of such spheroids 
today of p G = 10~ 3 Mpc~ 3 . 

Thus, luminous QSO MBHs at z ~ 3 are only about 10~ 3 as numerous as their galaxy 
host descendents today. One way to resolve this discrepancy is to assume that QSOs 
have very short duty cycles in their bright phase, of order t] ~ 10~ 3 . Since the quasar 
epoch runs about tq ~ 10 9 yr from z ~ 1.5 to z ~ 3.5 (the FWHM of the quasar plot 
in Figure 2), the lifetime in the bright phase would then be only t up ~ rjtg = 10 6 yr. 
The fractional mass change in this phase is only d~M/M = t up /ts ~ 0.02 eo.i, where 
the "Salpeter time" t$ (for an accreting black hole to e-fold in mass) is 

t s = M/M = 4 x 10 7 e .iyr, (5) 
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where we have parameterized the radiative efficiency in terms of eo.i = e/0.1 be- 
cause popular geometrically-thin optically-thick accretion disk models rarely exceed 
efficiencies of eo.i = 1. 

This leads to the disturbing conclusion that QSOs accrete only a tiny fraction of their 
mass while in their bright phase. This result, which depends critically on comparing 
the upper end of the AGN luminosity function with the upper end of the present- 
day MBH mass spectrum, is in sharp contrast with the near equality of the integral 
quantities (see Figure 1 & eqn 3): the total AGN background is roughly consistent 
with the total mass density in present-day MBHs, underpredicting the latter by a 
modest factor of 5 (consistent with significant accretion occurring as an advection 
dominated flow [[54]]). Since this integral constraint is dominated by bright objects, 
this conflict probably reflects a misidentification of the current fossil masses of the 
bright quasars. 

A plausible explanation is that MBHs may not have stayed constant in mass from 
the QSO era until now, but rather grew in mass by an average factor F. This growth 
might occur because hierarchical clustering merged these MBH's along with their 
protogalaxies. These mergers need not emit light. In that case, luminous quasars 
should be identified with MBHs today that are larger by the factor F, and their 
spheroids would be brighter by nearly the same factor. The exponential cutoff in the 
bright end of the luminosity function makes such spheroids much rarer and closer in 
abundance to the space density of QSOs. 

Specifically, for F = 5 the limiting MBH mass today becomes 2 x 10 9 M Q , the limiting 
spheroid luminosity is 4 x 1O 1O L , and the new benchmark galaxy in the luminosity 
function has Mb = —22.2, which is 300 times rarer. The difference in space density 
between quasars and spheroids is reduced to a factor of 3, the duty cycle comes up to 
1/3, and the lifetime of the bright phase is ~ 3 x 10 s yr. In this model bright quasars 
now spend a few Salpeter lifetimes in the bright phase, which is more plausible. 

Is growth by a factor F ~ 5 reasonable? Data on the recently discovered Lyman- 
Break Galaxies (LBGs) indicate slightly more growth. LBGs have been suggested by 
several authors as the early-formation phase of spheroidal components J55], [57[ . 
Small radii (1 — 2kpc |58], |56[| ) and small velocity dispersions (measured for only a 
handful of the brightest objects indicate modest masses of order a few xlO 9 M & . 
If these merge to form typical spheroids of today with L ~ 10 10 L & and masses 
M ~ 5 x 1O 1O M0, growth by more like a factor of 10 would be required. Models of 



hierarchically clustering protogalaxies [60] also suggest growth by F ~ 10. 



A second item that favors more growth is the lack of detection of AGNs in LBGs. If 
MBHs are forming everywhere in protogalaxies together with stars at the universal 
ratio M,/M star = 0.005 (Figure 1 ), then the average bright LBG should contain an 
MBH of about 1O 7 M , with Eddington luminosity Le ~ 10 45 ergs/sec. The apparent 
magnitude of such an object at z = 3 (including absorption by the factor 3-10 that 
affects LBGs ^TJ) is about tub = 26 — 27 , arguably faint enough to have escaped 
spectroscopic detection so far. 
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4.2 The Era of Galaxies 



If both MBHs and mergers are common among spheroids, two galaxies with pre- 
existing MBHs will frequently merge. The MBHs will spiral towards the center of the 
merger remnant, heating and perhaps ejecting stars from the center |T|. It is observed 
|62|1 that low-luminosity ellipticals and spiral bulges possess steeply rising central star 



profiles that approximate power laws, while high-luminosity ellipticals have central 
profiles that turn over to form much less dense centers (termed "cores"). The dense 
power-law centers of the small spheroids, their disky isophotes, and rapid rotation 
suggest that gas-rich, dissipative mergers or collapse may have controlled their mass 
distributions 0, [64], [65| , perhaps in the presence of one or more MBHs. 



The formation of the low-density cores of luminous core galaxies is more difficult. Cos- 
mological simulations tend to produce dark matter distributions with dense centers 
|66| , and dissipation is likely to sharpen their centers as would adiabatic compression 



by accretion of mass by a central black hole. pTJ, K5| In the dissipationless merger 



of gas-poor galaxies, the orbits of the associated MBHs would scour out a core in 



the stellar mass distribution, with a mass that tracks the black hole mass [67]. The 
predicted core masses and radii are a fair match to those observed, assuming that 
each spheroid has a BH given by the standard mass ratio M,/M staT = 0.005 ||65|| . 
The dissipationless mergers envisioned here are consistent with other characteristics 
of core galaxies such as their low rotation and boxy isophotes. An additional virtue 
of MBHs in core galaxies is that they may protect the cores from infilling by accreting 
low-mass, high-density satellites (by tidally shredding them), which may otherwise 
give luminous galaxies very bright centers. 



5 Tests of the Picture 

The detection of supermassive black holes and the discovery of dark matter share a 
common feature. In both cases there was skepticism of dynamical mass measurements, 
and acceptance was preceded by decades of debate. Dark matter has since come to be 
an essential feature of our understanding of many phenomena, ranging from galaxy 
rotation curves to the formation of structure in the Universe. We now seem engaged 
in a similar transition in the prevailing view of the centers of galaxies. 

An inevitable source of fuel for dead quasar engines is the debris from tidally dis- 
rupted stars. MBHs with masses M. < 1O 8 M disrupt main-sequence stars rather 
than swallowing them whole. Some of the debris from the star is ejected, but a por- 
tion remains bound to the MBH, forming an accretion disk that undergoes a "flare" 
lasting a few months to a year []2SL . Plausible models predict a V-band lu- 
minosity of about 1O 9 L . The event rate is controlled by how quickly stars can drift 
into the "loss cone" of low-angular-momentum orbits that come close to the MBH. 
Calculations show that faint, compact galaxies (e.g., M32) have the highest disrup- 
tion rates, about 10" 4 yr _1 . Larger, more diffuse galaxies (e.g., M87) have much lower 
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rates of about 10 -6 yr -1 , and often have sufficiently massive black holes to consume 
main sequence stars whole. It is possible that a stellar disruption by an MBH has 
already been witnessed spectroscopically. The nucleus of the spiral galaxy NGC 1097 
exhibited an ephemeral, broad, double-peaked Ha emission line [[nj, whose profile 
matched that expected from an accretion disk f7l |. 



Finally, it may be possible to detect the gravitational wave signature of merging 
MBHs, and thereby constrain the merger history of galaxies. In hierarchical models 
a typical bright galaxy has merged a few times since the quasar era. The timescales 
for decay of binary black holes in different regimes indicate that, for MBH with 
M.>10 7 M Q , the binary holes will merge on a timescale short compared to the next 
merger time |72"|. The merger rate for galaxies above 0.01L* may exceed 1/yr in the 



visible Universe. For an equal-mass binary black hole, the final orbit produces a lu- 
minosity of order L grav ~ c 5 /G = 10 60 ergs/sec in gravitational radiation, independent 
of M,. These mergers are the most powerful events in the Universe, but ironically 
they may not produce electromagnetic radiation. The energy is emitted over a time 
t oc GM./c 3 . The distinctive signature of a supermassive MBH merger as opposed to 
two stellar-mass black holes is the lower frequency and longer duration. Two 10 7 M Q 
black holes radiate much of their energy at a frequency about 10 -4 Hz. These events 
are too slow for LIGO (the Laser Interferometry Gravitational-wave Observatory, but 
are easy for LISA (the Laser Interferometric Space Array proposed as a Cornerstone 
Mission for ESA). A key test of the ideas in this paper is the observation of gravi- 
tational radiation from merging black holes at the centers of merged galaxies since 
z ~ 3. 
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